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ABSTRACT 

We report measurements of the millimeter wavelength brightness fluctuations produced by the 
atmosphere above the South Pole made with the Arcminute Cosmology Bolometer Array Receiver 
(Acbar). The data span the 2002 Austral winter during which Acbar was mounted on the Viper 
telescope at the South Pole. We recover the atmospheric signal in the presence of instrument noise 
by calculating the correlation between signals from distinct elements of the Acbar bolometer array. 
With this method, it is possible to measure atmospheric brightness fluctuations with high SNR even 
under the most stable atmospheric conditions. The observed atmospheric signal is characterized by the 
parameters of the Komolgorov- Taylor (KT) model, which are the amplitude and power law exponent 
describing the atmospheric power spectrum, and the two components of the wind angular velocity at 
the time of the observation. The KT model is typically a good description of the observed fluctuations, 
and fits to the data produce values of the Komolgorov exponent that are consistent with theoretical 
expectations. By combining the wind angular velocity results with measurements of the wind linear 
velocity, we find that the altitude of the observed atmospheric fluctuations is consistent with the 
distribution of water vapor determined from radiosonde data. For data corresponding to frequency 
passbands centered on 150, 219, and 274 GHz, we obtain median fluctuation power amplitudes of [10, 
38, 74] mK^ rad^^^^ in Rayleigh- Jeans temperature units. Comparing with previous work, we find 
that these median amplitudes are approximately an order of magnitude smaller than those found at 
the South Pole during the Austral summer and at least 30 times lower than found at the ALMA site 
in the Atacama desert. 

Subject headings: atmospheric effects — site testing — cosmic microwave background 



1. INTRODUCTION 

The geographic South Pole is one of the most im- 
portant sites for astronomical observations at mm and 
sub-mm wavelengths. This is primarily due to the rel- 
atively high altitude and the lack of water vapor in the 
atmosphere. In addition, the stability of the atmosphere 
makes it especially well suited for observations of the cos- 
mic microwave background (CMB) radiatio n. The recent 
detec tion of CMB polarization by DASI l)Kovac et al.l 
l2002|) an d fine scale CMB anisotropy measurements by 
Acbar ||Kuo et al.l 12004 hereafter K2004) precede a 
wave of new instruments to be deployed at the South 
Pole in the next few years. Compared to balloon-borne 
experiments, ground based experiments can make use of 
larger mirrors and achiev e longer con tinuous observation 
time. It was shown in iR.unvan ~a l. ( 2002J that the 
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background limited photon noise of the Acbar experi- 
ment is only 50% more than that a chieved in the ba lloon- 
bornc BOOMERanG experiment ijCrill et al.l l 2003D . For 
many types of observations at the South Pole, atmo- 
spheric brightness fluctuations will not significantly in- 
crease the noise above the background limit. 

The impact of atmospheric fluctuations on astronom- 
ical observations at mm wavelengths has been studied 
by several authors. Most of this work is concerned 
with distortions of the wavefront as the radiation prop- 
agates through the atmosphere and th e effects th is has 
on high resolution interferometry fe.g. lLavlll997|) . Pre- 
vious studies of atmospheric brightness fluctuations at 
mm and cm-wavelengths have been limited by sensitiv- 
ity ( WoUack et al. 1997; Lay & Halverson 2000, hereafter 
LH2000). Furthermore, this is the flrst characterization 
of the atmospheric brightness fluctuations made during 
the Austral winter at the South Pole. It is expected 
that the fluctuation power will be lower than during the 
Austral summer due to lower ambient temperatures and 
lower precipitable water vapor. 

We fit the observed atmospheric induced correla- 
tions between detectors in the AcBAR data with a 
Kolmogorov- Taylor (KT) model for the angular power 
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spectrum of atmospheric fluctuations. The high SNR 
with which we measure the atmospheric fluctuations 
makes it possible to determine the amphtude and ex- 
ponent of the KT power law and wind angular velocity 
for our atmospheric model. Combined with radiosonde 
measurements of wind linear velocity, we can infer the al- 
titude at which the fluctuation signals arise. This work 
provides a detailed characterization of fluctuations in 
the atmosphere above the South Pole that can be used 
to compare sites and predict the performance of instru- 
ments and observation strategies. 

In Section 12.11 we give a general overview of the en- 
vironment of the South Pole as pertains to millimeter 
wavelength astronomical observations. Section 12.21 de- 
scribes the spatial (Kolmogorov) and temporal (Taylor) 
components of the model. In Section |31 we detail the in- 
strument specifications and explain the observing strat- 
egy. Sections 0] and O discuss the methodology used to 
process the data and compare to models for the atmo- 
spheric fluctuations. In Sectional we present the results 
of this analysis, including the exponent of the power law, 
the amplitude of the fluctuations, measurements of the 
wind angular velocity, and the scale height of the fluc- 
tuations. Finally, in Section [7| we make some concluding 
remarks, including a comparison of our results with those 
obtained at the South Pole during the Austral summer 
and at the ALMA site in Chile. 

2. ATMOSPHERE ABOVE SOUTH POLE 
2.1. General Properties 

The large scale air flow at the South Pole is dominated 
by a katabatic (moving downward due to cooling) wind 
from the Antarctic plateau driven toward the coast. At 
the South Pole Station, radiosonde meteorology measure- 
ments are taken twice a day during the summer, and once 
a day in winter. These data, which include the pressure, 
temperature, relative humidity, and wind angular veloc- 
ity as a function of altitude, provide important informa- 
tion on the atmospheric properties. The weather data 
show that the surface wind direction is relatively con- 
stant, with higher wind speed in the winter than in the 
summer. It is also found that during the winter months, 
the cold surface creates a layer of temperature inversion 
where the air temperature rises from ~ — 60°C at ice 
level to ~ — 40°C at an altitude of approximately 200 m. 

Using balloon-borne microthermal probes, 
iMarks et alJ lIlQQQd) found this temperature inver- 
sion layer to be turbulent. The AASTO team also 
conflrmed this result with a non-intrusive SODAR 
system (Travouillon ct al. 2003). The lack of radiative 
heating during the winter and the featureless landscape 
near the South Pole imply that the turbulence in the 
inversion layer is due to shear-induced instability. The 
high correlation between the horizontal wind speed and 
temperature fluctuation reported by AASTO also sup- 
ports this point of view. The low-altitude atmospheric 
turbulence creates variations in the refracti on index, and 
is of primary concern for optical seeing l)Marktj f2002l: 
lTravouillon'eral] l200l . 

Observations at mm wavelengths are limited by the 
emission of radiation by molecules in the atmosphere. 
The primary components of the air (O2 and N2) are very 
well mixed; these "dry" components produce a uniform 



background of radiation that contributes photon noise 
to the detector. Water vapor in the atmosphere not only 
radiates at mm wavelengths, but also fluctuates in its 
mass fraction, resulting in brightness temperature fluc- 
tuations. Because of the low temperature immediately 
above the ground in the winter, the integrated water va- 
por pressure is small in the turbulent boundary layer 
(<200m). 

In the optically thin limit, the water vapor contribu- 
tion to sky brightness temperature is a function of the 
optical depth t and thermodynamic temperature T; in 
other words, it is proportional to the line of sight inte- 
gral of water vapor pressure. Approximating the satu- 
rated water vapor pressure as an exponential function 
in T (following the Goff-Gratch equation^), the water 
vapor pressure proflle can be derived from the publicly 
available archived radiosonde data.'^ Fig.nshows the wa- 
ter vapor profiles measured with radiosonde data during 
2002. Based on this data, we expect that the microwave 
radiation from water vapor at the South Pole comes pre- 
dominantly from a thick layer between 0.3 km and 2 km 
above the ice. 

One common way to characterize the quality of an ob- 
serving site is to report the distribution of the integrated 
precipitable water vapor for an extended period of time. 
Previous studies of sky opacity in the sub-mm produced 
median precipitabl e water vapor levels at the South Pole 
of about 0.25 mm l)Chamberlin fc BallvllTop) . The ra- 
diosonde data in Fig. ^ can be integrated to determine 
the precipitable water vapor for each day. Fig.Elshows a 
histogram of the precipitable water vapor measurements 
over the 2002 Austral winter. The mean precipitable 
water vapor during winter is 0.26 ± 0.2 mm. This is 
remarkably low when compared to the best six months 
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Fig. 1. — Water vapor pressure as a function of geopotential 
height above ice level at the South Pole as measured by radiosonde 
data, and normalized to the maximum observed value. Black lines 
represent water vapor pressure profiles for each day during the 
Austral winter; the white lino is the average of all the profiles. 
Most of the water vapor pressure originates in a broad layer in the 
atmosphere between altitudes of 0.3 and 2 km. Anomalous profiles 
well above or below the average appear to be due to measurement 
errors. 

http:/ /cires. colorado.edu/~voemel/vp. html 
^ ftp:/ /ice. ssec.wisc.edu/pub/raob/ 
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Fig. 2. — Distribution of integrated precipitable water vapor (in 
mm) at the South Pole for the 2002 Austral winter. The mean 
value of 0.26 mm is extremely small when compared to other well 
characterized astronomical sites. 



at other weh characterized sub-mm sites such as Mauna 
Kea (1.65 mm) and the Atacama desert (1.00 mm) l|Lanel 

mil). 

2.2. The Model 

The Kolmogorov Turbulence (KT) model describes the 
cascade of turbulent energy from the large scales on 
which it is input, to the small scales where it is dissi- 
pated. Large scale turbulence in the atmosphere above 
the South Pole is produced by processes such as wind 
shear and convection. These processes produce fluctua- 
tions with scales that exceed or are comparable to the 
thickness of the turbulent layer. Eventually, the turbu- 
lent energy is dissipated on scales comparable to 1 mm, 
which are far smaller than we are capable of observ- 
ing. The fluctuation angular power spectrum between 
the scales of energy injection and viscous dissipation, 
commonly known as the inertial subrange, is described 
by a power law with an exponent of —11/3. In the work 
described here, the scales being observed all lie in the in- 
ertial subrange and the fluctuation power should be well 
described by the KT model. 

We assume the atmospheric signal seen by the tele- 
scope is dominated by fluctuations in a layer of turbu- 
lence with thickness Ah, located at height h above the 
ground. The telescope is pointed at an elevation angle 
e, and the observed patch of sky subtends an angle of 
AO ^ 3°. Due to the relatively small size of the observed 
field and the typical altitude of the atmospheric water 
vapor. 

Ah > —AO , (1) 
sure 

and it is safe to assume that we are observing the pro- 
jection of the full 3-D Kolmogorov turbulence. In this 
limit, the thickness of the turbulent layer only affects 
the projected amplitude of the fluctuations. In addition, 
we are able to employ the small angle approximation 
{A6 <C 1) and avoid the complexity of spherical har- 
monic functions. Under these assumptions, the angular 
power spectrum follows the 3-D isotropic KT power law 



(P(fc) oc fc-"/3), and can be parametrized as 

P(a) = i3^in(6)-«/3(al + «2)-V2^ (2) 

where k is the spatial wave number, a = kh/ sine is the 
angular wave number, and is the amplitude fluctua- 
tion power at observing frequency v (in GHz) . From the 
KT model for atmospheric turbulence, we expect that 
b — 11/3; however, initially, we leave it as a free param- 
ete r in t he m odel fitting. This is the same model used 
bv ILavl l|1997f) and LH2000 to characterize atmospheric 
brightness fluctuations. is the amplitude of the power 
spectrum normalized to observations at the zenith and 
has units of mK^ rad"^^'^. We provide a more rigorous 
derivation of this relationship in Appendix ^ The (in- 
stantaneous) angular correlation function is the Fourier 
transform of the angular power spectrum P{a): 



C{9) = Jd^aP{a)( 



= / da dip aP{a)e 
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= 27r / da aP(Q;) Jo(27rtt 
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27Tia9 cos 



(3) 



To model the temporal variations, we fo llow the frozen 
turbulence hypothesis (FTH) proposed bv lTavloil l|l938|'l . 
which states that the small-scale turbulence velocity is 
much smaller than the global advection flow. In other 
words, the atmospheric fluctuations are modeled by a 
frozen pattern of fluctuations with Kolmogorov spec- 
trum, blown through the experimental beams by a uni- 
form wind. This model has been tested and verified by 
many groups engaged in interferometric observations. In 
the AcBAR observations, a chopping mirror sweeps the 
beam across the sky. In this case, applying FTH requires 
that the small scale turbulence velocity be much less than 
the sum of the chopper and wind angular speeds. 

We expect the large scale air flow over a flat terrain like 
the South Pole plateau to be nearly horizontal. Indeed, 
iTravouillon et al] (|2003*) report that the vertical wind 
speed is less than several tens of centimeter per second 
(compared to typical horizontal windspeeds on the order 
of 10 meters per second). In our model, we take the wind 
velocity v to be completely horizontal. 

3. OBSERVATIONS 

The data used in this study were collected using the Ar- 
cminute Cosmology Bolometric Array Receiver (Acbar, 
for a detailed description, see Runvan ct al. (2003)): an 
instrument designed to measure the anisotropy of the 
CMB. Acbar is an array of 16 bolometric detectors 
cooled to 240 mK with a 3 stage-'^He sorption refriger- 
ator. The square array is divided into four rows of four 
detectors each. The top and bottom rows correspond to 
274 GHz and 219 GHz (bandwidths 48 and 31 GHz, re- 
spectively), while the middle two rows are sensitive to 
150 GHz (bandwidth 31 GHz). Each detector is offset 
from its nearest neighbors by ~ 16'. AcBAR is mounted 
on the Viper telescope, a 2m off-axis Gregorian tele- 
scope at the geographical South Pole. The telescope 
produces beams with FWHM of 4 — 5' in all frequency 
bands. Through the use of a chopping flat mirror, the 
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telescope sweeps the beams across the sky at frequencies 
up to several Hz, with peak-to-peak amplitude of ~ 3° 
in the azimuthal direction. The results presented in this 
paper are derived from the data taken with a chopping 
frequency of 0.3 Hz. Acbar has high resolution in spher- 
ical harmonic multipole Z-space {Al ^ 100) over a wide 
range of angular scales (150 < ^ < 3000), allowing precise 
measurements of the CMB power spectrum (K2004). 

The observations used in this work began in late March 
of 2002 and continued through the end of July when 
the supply of liquid helium at the South Pole station 
was exhausted. There are about 200 data files spanning 
that period, each of which constitutes up to six hours of 
continuous observation. One six hour observation con- 
sists of approximately 300 "stares" each lasting about 
one minute. During a stare, the telescope tracks a fixed 
position on the sky while the chopper sweeps out a trian- 
gle wave in azimuth. Each of these stares falls into one of 
three categories: lead, main, or trail (LMT). Lead and 
trail stares are offset ±3° in RA from the main stares 
and have half the integration time that main stares do. 
Each LMT group of stares is offset in declination, so that 
the fin al dataset represents a rastered map of the sky. 
iKuo et al. ( 2004) produced a CMB angular power spec- 
trum from the M-(L+T)/2 differenced map. However, 
for the purpose of this atmospheric study, LMT stares 
are not differenced and are treated identically. 

4. DATA PROCESSING 

The beams of the Acbar array are swept in azimuth 
with a 0.3 Hz triangular oscillation. The amplitude of the 
waveform is 3° in azimuth when the telescope is pointed 
at the horizon. The timestream from the Acbar de- 
tectors is binned into vectors, or data strips, , i = 
1,2, ...,n, each of which consists of n = 128 tempera- 
ture measurements for array element p covering 3° on 
the sky. The data corresponding to sweeps in the right 
(positive AZ) or left (negative AZ) directions are binned 
separately, resulting in two data vectors from every com- 
plete chopper cycle. The data vectors are calibrated in 
RJ temperature units by referencing them to observa- 
tions of the br ightness temperature of Mars and Venus 
ijRunvan et alJi2003.) . 

Correlation matrices are calculated from the outer 
products of data vectors corresponding to different el- 
ements in the detector array and compared with the 
model for the atmosphere-induced correlations. By cor- 
relating data vectors from distinct array elements, we ob- 
tain cross-correlation matrices that have a significantly 
lower bias from noise power than the autocorrelation. To 
maximize the signal-to-noise ratio (SNR), we correlate 
all pairs of array elements within a particular frequency 
passband. The angular offsets between these pairs of ar- 
ray elements are listed in Tabled and a schematic dia- 
gram of the pairs is presented in Fig.|31 Because there are 
twice as many array elements operating in the 150 GHz 
band, there are many more correlation matrices for that 
frequency passband. 

The statistical correlation depends only on the relative 
displacement of the pairs of array elements. Therefore, 
we average the correlation matrices of all pairs with the 
same frequency band and displacement vector together 
in a single correlation matrix. In order to reduce the 
volume of data, we assume that the correlation is sta- 
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Note. — Angular offset in azimuth (A(^) and elevation (Ae) be- 
tween array elements used for the correlation analysis. The total 
number of correlation maps associated with each array displace- 
ment and chopper direction is listed for the 150 GHz (two rows of 
detectors) and 219, 274 GHz (one row of detectors each) data. The 
219 and 274 GHz data contain only displacements in azimuth and 
have the same number of correlations maps; 71219,274 corresponds 
to either of the two detector rows. 



tionary and average the matrices over a period of many 
sweeps. The rotation of the telescope with respect to 
the ground based wind vector places constraints on the 
acceptable period of time averaging. The fundamental 
length of time over which the correlations are averaged 
is one quarter of a data file or about 1.5 hours (see Sec- 
tion for details). This procedure is repeated with 
the data from the opposite direction of chopper motion 
to yield two correlation matrices for each band, detector 
displacement vector, and time period. After removing 
modes that are susceptible to airmass gradient and chop- 
per synchronous offsets (zeroeth, first, and second order 
polynomial) and subtracting the average of all the stares 
in the file from each stare, we compare the measured 
correlation matrices with those calculated from eq. jS] . 

The chopping mirror on the Viper telescope sweeps the 
optical beams at a constant angular speed ±7 (rad s~^). 
At time zero (t = 0), a detector element measures sky 
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Fig. 3. — Schematic of ACBAR focal plane demonstrating an 
example of each array element pair correlation that has been used 
in the analysis. See Table for a tabular representation of all 
possible pairs. Correlations are computed for all possible pairs of 
distinct array elements within a frequency passband. The letters 
A, B, C ... correspond to the array displacements in Table HI 
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temperature Ti in the direction (AZ, EL)= (cj) ,£). r sec- 
onds later, another detector element (with position offset 
A(j), Ae) measures sky temperature T2 in the direction 
of ((/)', e'). The correlation of Ti and T2 is C{6) given by 
eq. P], where the angular separation 9 is 



9^^{^' -cby + ie' -ey (4) 

= \J {w^T TIT - Acjyy + [w^T - Ae)2 , 

and {w^,We) are the components of the wind angular 
velocity in the direction of (0, e). The wind linear veloc- 
ity is related to the angular velocity by: 



h 



Sine 



W0 ■ 



sm e 



sin^ e 



(5) 



Note that the angular separation 9 is measured on the 
"frozen" screen of atmospheric fluctuations. When A(j> — 
Ae = 0, C{9) becomes the autocorrelation function. One 
interesting implication of this scan strategy is that the 
observed atmospheric fluctuation signal depends on the 
direction of the chopper motion (a plus or minus sign in 
eq. ^) with respect to the wind direction. Compared to 
a drift scan experiment (7 = 0), the fluctuation signal 
is suppressed when the beams and the wind move in the 
same direction, and enhanced when they move in the op- 
posite directions. This effect is easily understood in the 
limit where 7 = w^, We = 0. The optical beams follow 
the frozen turbulence as it is blown across the sky and the 
AC-coupled bolometers will not detect any atmospheric 
signal. On the other hand, when the chopper goes in 
the opposite direction to the wind, the instrument will 
record an enhancement in atmospheric signal. 

5. MODEL FITTING 

The KT model given by eq. [2] is used to character- 
ize the atmospheric fluctuation angular power spectrum. 
For each observation, we determine a set of parameters 
that minimizes the difference between the theoretical cor- 
relation function and the observed data. These parame- 
ters are defined in eqs. [21-E] and include the amplitude 
for each observing frequency, B^; a power law exponent, 
6; and angular wind velocity in the x and y directions, 
Vx/h and Vy/h, respectively. As a result of the Earth's 
rotation, the telescope tracks across the sky in azimuth, 
changing by 90° during the course of a six hour observa- 
tion. Using eq. [S], we find the components of the pre- 
vailing wind angular velocity, v^/h and Vy/h as functions 
of the wind velocity components in the telescope frame 
and the time during the observation (which sets 0, the 
azimuth) . 

To compare the parameterized model to the data, we 
first need to generate a theoretical correlation matrix, C, 
from eq. P] . We do this by discretizing the temporal lag 
T in eq. 0] into n bins, where n is the length of each data 
vector being correlated — in our case, 128. If we denote 
Ti as the discretized temporal lag associated with the sky 
temperature measurement in the direction of (0 , e) and 
Tj as the same quantity for a distinct detector in the 
direction of (0', e'), then the temporal lag between these 
two temperature measurements is r = — t, and we can 



generate the discrete, n x n correlation matrix C from 
the correlation function as follows: 



C'^ = C(6I(t, A0,Ae)). 



(6) 



where Aip = (/>' — (p and Ae — e' — e. We perform an 
equivalent mode removal on C*^ as was performed on the 
observed data, so that the theoretical and observed cor- 
relations can be directly compared. We calculate a 
for the fit of the theoretical correlation matrix to the 
observed data for each frequency passband and unique 
correlation matrix: 



Xk,. 



'del 



data 



(7) 



where ak,ii represents the intrinsic detector noise asso- 
ciated with the k*^ correlation map for observation fre- 
quency ly. We compute the correlation matrices corre- 
sponding to the noise in a given 6 hour observation (de- 
noted by Ck^v^noiae) described in Section ^ with the 
difference that the data vectors for the detector pairs are 
drawn from the three possible pairs of 1.5 hour subfiles 
separated by at least 1.5 hours in time. We expect the 
atmosphere to be uncorrelated between these two mea- 
surements and the correlation should provide an accurate 
estimate of the noise. Finally, ak,ii is obtained from the 
variance of the /c*'' correlation matrix of frequency i^: 



(8) 



The overall used to determine the model parameters 
is given by 



(9) 



where Wk^v represents the weight of the fc*'' correlation 
matrix found in column A {v = 150 GHz) or 5 (j^ = 219 or 
274 GHz) of Table H] For example. A: = 0, = 274 GHz 
corresponds to left- or right-going correlation matrices 
with angular offsets of A0 = —16' and Ae = 0'. From 
Table ^ we see there are three pairs with this angular 
offset, meaning wo,274 = 3. k — X^v — 150 GHz corre- 
sponds to left- or right-going correlation matrices with 
angular offsets of A0 — —32' and Ae = 0', so Wi.i50 = 4. 

The lack of detector differences in the e direction for the 
219 and 274 GHz data results in a poor determination of 
this component of the wind velocity from this data alone. 
However, the 150 GHz data feature correlation matrices 
with angular offsets in both (j) and e directions, and can 
be used to obtain reliable wind velocity parameters. We 
simultaneously fit the data for all frequencies with the 
wind velocity to be constrained to be the same for all 
detector pairs. The 150 GHz data provides constraints 
on the e component of the wind velocity. We use the 
IDL AMOEBA routine to perform the minimization and 
determine the best fit parameters for the KT model. 

5.1. Correlation Matrix Fits 

We begin by fitting the model parameters to the data 
divided in 1.5 hour segments. Wc find that the atmo- 
spheric fluctuations observed by Acbar are well de- 
scribed by eqs. [SJ 0] for about 93% of the data. A 
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Fig. 4. — Correlation matrices produced by the model (1st and 3rd columns) as compared to the observed correlation matrices (2nd and 
4th columns) for a 1.5 hour file of 150 GHz data from May 30, 2002. The two columns on the left correspond to a left-going motion of the 
chopper while the two columns on the right correspond to a right going sweep. The rows are a representative sampling of available angular 
displacements between 150 GHz array elements. 



typical example of a good fit by the model is shown in 
Figs. 21 and [S] These figures show the theoretical and 
observed correlation matrices for a f .5 hom' data file ob- 
tained on 2002 May 30. The two columns on the left 
correspond to left-going sweeps of the telescope, while 
the two columns on the right represent the opposite di- 
rection. Fig. [51 shows a cross-section through the center 
of the correlation matrices presented in Fig. 0] This file 
has a particularly high azimuthal component of wind ve- 
locity and the cross-section plots demonstrate the differ- 
ence in amplitude between the left-going and right-going 
sweeps. As described in Section 0] the sky noise is sup- 
pressed when the beam follows the wind (as seen in the 
left-going sweeps for this day), and enhanced when they 
move in the opposite directions. 

One of the unique features of this work is that it ap- 
pears possible — provided the KT model reasonably de- 
scribes the data — to measure the atmospheric fluctua- 
tion power under nearly all atmospheric conditions, even 
during the best weather. As mentioned in Section [S] 
the intrinsic detector noise sets a fundamental limit to 
the smallest atmospheric fluctuation power that can be 



measured. We characterize the noise in the fluctuation 
amplitude by fitting to the random correlations between 
detectors in time-separated sweeps. In Fig.|Sl we plot the 
noise power amplitude vs. the power amplitudes found 
for the correlations within the sweeps. The line repre- 
sents a SNR of approximately one. Nearly all of the 
data lie above this line, indicating that the amplitude of 
the fluctuations is not limited by an instrumental noise 
floor and that the amplitude of the fluctuation power is 
typically measured with high accuracy (median SNR is 
140) under nearly all observing conditions. 

We identify the observations that do not provide ade- 
quate fits to the KT model from the defined by eq. [5] . 
We find that 7% of the total observations fail to measure 
the KT model amplitude with < 20 * median{x^)- In 
general, the poor fits are not due to a lack of fluctuation 
power, but rather the presence of fluctuations that can- 
not be described by our assumed model. We have broken 
these files into smaller time periods and find no improve- 
ment in the fits, indicating that it is not a problem with 
the fit parameters changing in time. Two possible expla- 
nations are that either the observed spatial fluctuations 
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Fig. 5. — Cross-Sectional plots of the observed correlation matrices shown in Fig. I^with the model fit overplotted. The cross-section is 
taken in the horizontal direction through the center of the correlation matrix, and serves as a demonstration of the quality of the fit. Due 
to the high wind velocity, the sky noise is suppressed for the chopper moving in the same direction as the wind (left column) or enhanced 
for motion in the opposite direction (right column). 



do not obey a Kolmogorov power law at this time, or the 
atmospheric structure is not frozen on the timescale of 
a sweep. Another more likely explanation is that a sig- 
nificant fraction of the observed power is coming from a 
layer in the near field of the telescope, which we estimate 
to be approximately 500 m. In this case, the assumptions 
made in the derivation of the theoretical correlation func- 
tion are not valid. The data with poor fits are typically 
associated with very large fluctuation amplitudes, and 
are therefore likely to be due to a severe storm. 

5.2. Timescale of Disturbances 

By averaging the correlation matrices over longer pe- 
riods of time, it is possible to measure the model pa- 



rameters (particularly the wind velocity) with higher ac- 
curacy. However, it is important that the atmospheric 
conditions remain constant over the period for which the 
correlations are averaged. The CMB observations we are 
using for this analysis cover 6 hours; this sets a practi- 
cal upper limit to the averaging. In this section, we show 
that the atmospheric conditions are typically quite stable 
over the course of a six hour observation. 

Each of the six hour data files is subdivided into four 
1.5 hour subfiles. Model parameters are determined for 
each of the subfiles from fits to the 150 GHz data. To 
reduce computing time, we fit only to the 150 GHz data 
for the tests described here and in Section [?31 In Fig.[7| 
we show the correlation between the power amplitudes 
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Fig. 6. — A graphical representation of the SNR in the de- 
termination of the fluctuation power amplitude. On the y-axis, 
we plot the mean fluctuation amplitude found from the individual 
1.5- hour segments of a 6 hour data file. The x-axis is the ampli- 
tude computed from the correlation between detector timestreams 
separated by at least 1.5 hours within a given 6-hour observation. 
The diagonal line represents equal fluctuation power amplitude for 
the signal and noise correlations. The fluctuation amplitude is re- 
covered with high SNR in nearly all files. In fact, the median SNR 
is 140. 



determined for the 1st and 4th quarters of each file. The 
strong correlation indicates that the atmospheric fluctu- 
ation power is relatively consistent in time and we are 
justified in fitting a single fluctuation power to a six 
hour data set. Due to the rotation of the telescope with 
respect to the prevailing wind direction as it tracks a 
source, we cannot simply average the correlation matrices 
over the entire 6 hour observation. Instead, we project a 
constant ground-based wind velocity into the appropriate 
telescope frame for each of the quarters of the flle. The 
minimization for the model parameters then proceeds as 
for a single flle; however, we use the combined for all 
four flies as a metric of goodness of fit. 




Fig. 7. — Correlation of brightness fluctuation power amplitude 
over the span of a six hour observation. We show correlation be- 
tween ^150 measured from first and fourth quarters of a six hour 
data file. The strong correlation implies that that brightness fluc- 
tuation power typically remains constant over a time period of at 
least six hours. 



5.3. Power Law Exponent 

We can determine the exponent of the KT power law 
for those observations in which we measure the fluctua- 
tions with high SNR. In Fig. El we show the distribution 
of measured values of the power law exponent over the 
course of the winter where the mean exponent for that 
period of time is & = 4.1 ± 0.8. If we restrict the observa- 
tions to those with above-median SNR, we find that the 
distribution of exponents is described by 6 = 3.9 ± 0.6. 
In both cases, our results are consistent with those ex- 
pected for the KT model in the 3D regime. For the rest of 
our analysis, we have set the exponent of the KT power 
law equal to 11/3. This assumption is physically moti- 
vated and increases the stability of the model-fitting by 
removing one of the free parameters. 

6. RESULTS 

In this Section, we describe the basic results of the fits 
of the KT model to the observed data. As described 
in Section 15.21 the stability of the atmospheric condi- 
tions allows us to determine the model parameters for 
an entire six hour observation. Although the fluctua- 
tion amplitude is determined with high SNR in each of 
the flle quarters, fltting to all four simultaneously signif- 
icantly improves the determination of the KT exponent 
and wind velocity components. Unless otherwise noted, 
the following results are found assuming that the model 
parameters are constant over a six hour data flle. The 
amplitudes for the fluctuation power are free to vary in 
each of the bands, while the wind angular velocity is con- 
strained to be constant across all bands. For the 219 and 
274 GHz rows, this greatly improves the stability of the 
results. 
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Fig. 8. — Distribution of the measured KT power law exponent 
values for the Austral winter of 2002. Observations that do n ot fit 
the KT power law (according to the cut described in Section [5.1i 
are not included in this distribution. The solid line shows the 
distribution of exponents for all good fits to the KT model; the 
mean of this distribution is 4.1±0.8. The dotted lines shows the 
distribution of exponents for all observations with below-mean 
for the fit to the KT model; the mean of this distribution is 3.9it0.6. 
In both cases, the results agree with the theoretical KT power law 
for small angular scales, which is shown as a vertical line. For the 
rest of the analysis, we fix the exponent at the expected value of 
b = 11/3 to remove the degeneracy between the amplitude and 
exponent of the power law. 



South Pole Atmosphere 



9 



6.1. Amplitude Measurements 

We have used fits to the KT model to measure the 
brightness fluctuation power amplitude of the atmo- 
sphere over a four month period from late March through 
July 2002. Although the measurements were made at a 
range of elevations, the results have been normalized in 
terms of observations at the zenith. The results of this 
analysis for 150 GHz data are shown in Fig. ^ The 
atmospheric brightness fluctuation power is relatively 
stable, with roughly 60% of the data falling between 
Bf;^Q = 5 - 40mK^ rad"^/^. This is in contrast to the 
findings of LH2000, who reported a bimodal distribution 
of brightness fluctuation power for the Austral summer 
of 1996. Solar heating and associated weather during the 
summer is one possible cause of this difference. 

We have performed a similar analysis for the 219 and 
274 GHz data. The emissivity of water vapor increases 
as a function of frequency, resulting in higher amplitude 
fluctuations at higher frequencies. In Fig. 1101 we com- 
pare the amplitudes found from the 219 GHz observa- 
tions with those found at 150 GHz. The degree of corre- 
lation is a consequence of the high SNR of all the mea- 
surements and is similar to that found between the 150 
and 274 GHz data. The brightness fluctuation power am- 
plitude at 274 (219) GHz is approximately 7.4 (3.7) times 
greater than at 150 GHz. In Fig. ^] we show a Cu- 
mulative Distribution Function (CDF) of the brightness 
fluctuation power for each of the three frequency bands. 
From these distributions, we extract upper limits to the 
quartiles in brightness fluctuation power amplitude and 
list them in Table |21 We compensate for the 7% of the 
data that do not fit the KT model by assuming that all 
the poor KT fits have amplitudes larger than the largest 
good KT fit. The 27th, 54th, and 81st percentiles then 
serve as upper limits to the quartiles of the full distribu- 
tion. 
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Fig. 10. — The fluctuation power amplitudes determined at 
150 GHz vs. 219 GHz. The high degree of correlation is a result of 
the large SNR in the determination of the amplitudes and the ro- 
bust fits to the model. A similar plot for the 150 GHz vs. 274 GHz 
results is qualitatively identical. 



The median fluctuation power amplitudes for each of 
the bands are listed in Table |31 For comparison, we 
also list the median fluctuation power measured by the 
Python experiment from the paper of LH2000. Due to 
the increase in atmospheric brightness with frequency 
and the higher sensitivity of the instrument, the ACBAR 
measurements are signiflcantly more sensitive to water 
vapor than those made with Python. This sensitivity 
is helpful in characterizing the much more quiet winter 
atmosphere. 

We have used a model for the spectrum of atmospheric 
water vapor to compare the results from the four fre- 
quency bands covered by Acbar and Python. The At- 
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Fig. 9. — The brightness fluctuation power, in mK^ rad"^/"^, for 
4 months of 150 GHz ACBAR observations at the South Pole. Each 
data point represents about six hours of continuous observations 
normalized to observations at the zenith angle. Approximately 7% 
of the data exhibit anomalous turbulence which are not well fit by 
a KT model, and are not shown here. We assume that these data 
always have amplitudes larger than data which can be modeled 
with a KT power law and adjust our amplitude distributions ac- 
cordingly. Mechanical problems with the instrumentation limited 
the acquisition of data for days 100-112. 




B/ [mK' rad-=''"] 



Fig. 11. — The cumulative distribution functions of brightness 
fluctuation power amplitude (in Rayleigh-Jeans temperature units) 
for the entire winter season of ACBAR observations. The solid, dot- 
ted, and dashed lines represent 150 GHz, 219 GHz, and 274 GHz 
data, respectively. The 7% of the data that is poorly fit by a KT 
model is not included in this plot. The upper limits to amplitude 
quartiles given in Table 1^ have been calculated from these distri- 
butions. 
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mospheric Transmission (AT)^ code is used to produce 
spectra corresponding to winter conditions at the South 
Pole with 0.25 mm water vapor as weU as 0.0 mm water 
vapor. The difference of these two spectra gives an esti- 
mate of the emission due to water vapor as a function of 
frequency. The emission spectrum of atmospheric water 
vapor is shown in Fig. ^1 with the four passbands shown 
in gray. We integrate over the passbands to obtain the ef- 
fective opacity r^aj^j due to 0.25 mm water vapor for each 
of the observing frequencies; 40, 150, 219, and 274 GHz. 
We then scale the fluctuation power amplitude for each 
of the bands by the square of the ratio of the water vapor 
opacities to normalize the results to the equivalent power 
that would have been seen in simultaneous observations 
at 150 GHz. 

Equivalent Bf^Q = (^^^ 

Given a perfect understanding of the Acbar bands, an 
accurate water vapor spectrum, and the validity of the 
assumption that the atmosphere is optically thin, we ex- 
pect the scaled fluctuation amplitudes for the Acbar 
measurements to be identical. The equivalent bright- 
ness fluctuation power amplitudes for the Python obser- 
vations made during the Austral summer are more than 
10 times higher than those for the Acbar measurements 
made during the Austral winter. We attribute this dif- 
ference largely to the lower water vapor content of the 
atmosphere during the winter, although the stability of 
the atmosphere may also be improved in the winter. 

The measurements of the fluctuation amplitude made 
with the Python experiment and reported in LH2000 are 
given in terms of Ah^^^, while the Acbar results are 
given in terms of B^ — Ah^^^. In order to compare the 
results of the experiments, it is necessary to divide the 
numbers from LH2000 by the altitude of the layer con- 
taining the fluctuations. This additional factor of h was 
introduced by the authors of LH2000 so that the data 
would be easier to insert into a set of approximate an- 
alytic expressions predicting the level of residual atmo- 
spheric noise in different classes of observations. How- 
ever, these expressions are only valid in the limit of long 
averaging times, which are not appropriate for many ob- 
serving strategies. B^ is what is actually measured by 
the experiments, and is therefore more fundamental and 
useful than Ah^^^. In Section IF?^ we show that the scale 
height of the fluctuations is consistent with the distribu- 
tion of precipitable water vapor. Should one desire to 
compute Ah^/^ from the Acbar data, multiply B^ by 
the water vapor pressure weighted altitude, h — 1300 m. 

Measurements of the sub-mm atmospheric opacity are 
an important monitor of observing conditions. We have 
investigated the correlation of the measured fluctua- 
tion amplitude with both the sub-mm opacity (r) and 
the variance of the opacity (cr,-). The sub-mm tipper 
sky brig htness monitor) on the AST/RO experiment 
Radfordll200^ at the South Pole measures the 350 /im 
zenith opacity and atmospheric temperature every 15 
minutes. The 350 /im atmospheric opacity is much larger 
than 150 GHz opacity, but they are both dominated by 
water vapor and, therefore, scale nearly linearly. We bin 

® Airhead Software, Boulder CO 



TABLE 2 



Percentile 


27th 


54th 


81th 


150 GHz 


3.7 


10. 


37. 


219 GHz 


11. 


38. 


160 


274 GHz 


28. 


74. 


230 



Note. — Upper limits to quartiles of fluctuation amplitude 
(BI) for 150, 219, and 274 GHz data in Rayleigh- Jeans temper- 
ature units of mK^ rad"^"^^. These values are calculated from 
the cumulative distribution functions shown in Fig. lllj We report 
the 27th, 54th, and 81st percentiles, because the distribution from 
which they are calculated does not include the 7% of non-KT ob- 
servations. We assume that each of the non-KT observations has 
an amplitude larger than the largest KT observation, the values 
in this table represent upper limits to the 25th, 50th, and 75th 
percentiles of the full distribution. 



the opacity values from the tipper to match the length of 
the observation (typically about six hours); the fractional 
variance in opacity over that time-span (ctt/t) is usually 
less than 10%. Fig. 1131 shows the fluctuation power am- 
plitude versus the square of the opacity, r^. The corre- 
lation between the fluctuation power amplitude and the 
variance of the opacityis similarly weak. These results 
imply that r and CTt are not, by themselves, an adequate 
judge of the stability of the atmosphere on timescales 
relevant to the Acbar measurements. 

6.2. Wind Velocity and Scale Height 

In general, the angular speed of the chopper is higher 
than the effective angular speed of the fluctuations pro- 
duced by wind. Therefore, the effect of the wind velocity 
on the observed correlations is typically sub-dominant 
to the effect of the chopper motion. In practice, the am- 
plitude of the fluctuations can usually be recovered with 
high SNR even when very poor constraints can be placed 
on the wind velocity. However, in cases where the fluc- 
tuations are measured with high SNR, the wind angular 
velocity can be accurately recovered. This wind angular 
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Fig. 12. — Water vapor emission spectrum obtained from the 
difference between a "wet" spectrum with 0.25 mm input water 
vapor and a "dry" spectrum with 0.0 mm water vapor. Approx- 
imate frequency passbands of the Python (40 GHz) and ACBAR, 
(150, 219, 274 GHz) experiments are shown with grey rectangular 
bars. 
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TABLE 3 



Frequency 40 GHz 150 GHz 219 GHz 274 GHz 

Median [niK^ rad-5/3] (CMB units) 0.90 31. 350 2.0x10^ 

Median [niK^ rad-^/^j (RJ units) 0.84 10. 38. 74. 

Equivalent B3 for 150 GHz 11. 1.0 0.92 0.86 



Note. — Comparison of median brightness fluctuation power for measurements with the AcBAR 
and Python (see LH2000) experiments. The results are listed in both Rayleigh- Jeans and CMB tem- 
perature units. Clearly, RJ units are more natural for comparing between different bands, however, 
we include the CMB normalization in order to facilitate the prediction of noise in future exper- 
iments. We use a water vapor spectrum to normalize the atmospheric fluctuation power for the 
Python 40 GHz and AcBAR 219, 274 GHz observations to the equivalent power for the 150 GHz AcBAR 
band. The data for the three AcBAR bands were taken simultaneously and, therefore, the equivalent 
150 GHz amplitudes should be roughly identical. The Python numbers are computed from Table 1 
in LH2000 assuming the average altitude of the fluctuations, /lat; — 1000m. The order of magnitude 
higher equivalent fluctuation power for the extrapolated 40 GHz Python results is a result of the lower 
precipitable water vapor and increased stability of the atmosphere during the winter at the South 
Pole. 



velocity can be combined with radiosonde measurements 
of wind linear velocity to estimate the scale height from 
which the atmospheric fluctuation power originates. 

Data files with a (as described by eq. below the 
median are selected as being likely to result in reliable 
wind velocities. The wind angular velocities for each 1.5- 
hour file are found in a coordinate system that is fixed 
with respect to the direction the telescope is pointing. 
However, over the course of a long data file, the telescope 
tracks a source through a considerable range of angle and 
this motion must be taken into account. The transfor- 
mation described by eq. |3] ties the local wind velocities 
for each 1.5-hour data set into a single wind velocity with 
respect to the ground constant over the entire data file. 
We are justified in doing this because both the measure- 
ments of the wind angular velocities with Acbar and 
the wind linear velocities with radiosondes show that the 
wind is typically quite constant over a 6-hour data file. 

The Acbar measurements of wind angular velocity 
can be combined with radiosonde linear velocity and wa- 
ter vapor pressure as a function of elevation to determine 
the mean elevation (h) from which the fluctuations orig- 



inate: 

h = (10) 

^Acbar 

where Vrs is defined as 

_ v{z)P^uv{z)dz 

Jo Pwv{z)dz 

v{z) is the wind linear speed from the radiosonde data, 
Pwvj z) t he water vapor pressure (calculated in Sec- 
tion [OJ, and z the altitude above ice level. Fig. 
is a histogram of the altitude of the inferred layer giv- 
ing rise to the measured fluctuations. The results are in 
reasonable qualitative agreement with the average wa- 
ter vapor pressure weighted altitude calculated from ra- 
diosonde data shown in Fig. ^ This provides further 
confidence that we understand the source of the fluctua- 
tions and have accurately characterized them. 
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Fig. 13. — Fluctuation power amplitude versus 350 /jm opacity- 
squared. The weak correlation between SfsQ and r'^ implies that 
T alone is not a reliable indicator of atmospheric fluctuations. 
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Fig. 14. — Distribution of the inferred altitudes for the atmo- 
spheric fluctuations. These results are computed from the ratio of 
the Acbar (below-median 150 GHz data) wind angular veloci- 
ties to the radiosonde derived water vapor pressure weighted wind 
linear velocities. The results are consistent with the average water 
vapor pressure profile calculated from the data shown in Fig.lTland 
reproduced here as a smooth line. The mean inferred scale height 
of the emission is ^ 1.3 km. 
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7. CONCLUSIONS 

We have produced a detailed characterization of the 
millimeter wavelength atmospheric fluctuations at the 
South Pole during the Austral winter of 2002. By corre- 
lating the signals between pairs of elements in the Acbar 
bolometer array it is possible to reduce the effect of 
uncorrelated detector noise such that the atmospheric 
fluctuations can be measured with high SNR even for 
the lowest fluctuation amplitudes. We found that a KT 
model is an accurate description of the atmospheric fluc- 
tuation power ~ 93% of the time. The median expo- 
nent of the power law is found to be consistent with 
the theoretically expected value b = 11/3. For the bulk 
of the analysis, we fix the exponent to 11/3 in order 
to improve the SNR in the determination of the other 
model parameters. We found upper limits to the the 
quartiles of the 150 GHz atmospheric fluctuation power 
for the four winter months of observation to be [3.7, 10., 
37] niK^ rad~^/^. The scaling between the measured fluc- 
tuation power amplitudes for each of the three AcBAR 
frequencies is consistent with that expected for emission 
from water vapor. We used this fact to compare our 
results with those of LH2000 made at 40 GHz with the 
Python experiment in the Austral summer of 1996. As- 
suming that the scale heights of the fluctuations in the 
summer and winter are equivalent, the median fluctua- 
tion power found with Acbar observations in the Aus- 
tral winter corresponds to a fluctuation power amplitude 
~ 10 times smaller than that found with Python during 
the Austral summer. From the analysis in LH2000, this 
implies that the atmospheric fluctuation power at the 
South Pole in the Austral winter is at least ~ 30 times 



smaller than that at the ALMA site in the Chilean An- 
des. In addition, we found that the fluctuation power am- 
plitude is only weakly correlated with atmospheric opac- 
ity or its variance on 15 minute timescales. Therefore, for 
experiments sensitive to atmospheric fluctuations, peri- 
odic measurements of optical depth do not appear to be 
an adequate characterization of an astronomical site. Fi- 
nally, we determined the wind angular velocity for each 
of the observations with below median and compared 
our value to the wind linear velocities obtained by ra- 
diosonde measurements to determine the scale height of 
the observed atmospheric fluctuations. We find that the 
inferred scale heights for each of the data files roughly 
fall within the 300 m-2000 m range suggested by the ra- 
diosonde water vapor pressure profiles. 

We believe this work provides a fairly complete pic- 
ture of the fluctuation power in the atmosphere above 
the South Pole. Already these results are being used to 
compare potential observation strategies for the 10-meter 
South Pole Telescope, which will be deployed in 2007. We 
anticipate that similar studies to that presented here, will 
be performed with data from other array receivers in or- 
der to compare quantitatively the atmospheric stability 
of observing sites. 
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ment. We would also like to thank Nils Halverson and 
Tom Crawford for their careful reading of an early draft 
of this work and their helpful suggestions for it's improve- 
ment. This work has been supported by the NSF office 
of polar programs grant OPP-0091840. 



APPENDIX 
PROJECTED POWER SPECTRUM 

Here we derive the instantaneous angular power spectrum, given the three dimensional temperature-density fluctua- 
tion spectrum. The thermodynamic temperature of the atmosphere is on the order of 200 K, hence the Ray leigh- Jeans 
law is a good approximation at mm wavelengths. In the optically thin limit and under the condition of local thermody- 
namic equilibrium, the speciflc radiation is proportional to the line of sight integration of the light emitting elements. 
Therefore, we can define the effective temperature TE(r), such that the brightness temperature Tb along f = r/r is 
given by 

Tb = J Tek dr , 

where k is the opacity/thickness of the turbulent layer for millimeter wave photons. We assume that Te includes 
both the physical temperature variations and density variations. The angular correlation function for Tb can then be 
written as 

<.(/i+A/i)/sine (/i+A/i)/ sin e 

C{d) = {TB{e) ■ Tb{0)) = dz dz'K^TE{x,y,z')TE{Q,0,z)) 

J h/ sin € J h/ sin e 

{h+Ah)/ sine ^{h+Ah)/ sine 

dz / dz' K^C3D{x,y, z' — z) . 

h I sin e J h/ sin e 

Fig. lAlSl shows the coordinate system used in this calculation. In the small angle regime, the angular power spectrum 
is the Fourier transform of the angular correlation function, which can be derived from the 3-D correlation function, 
which in turn is the 3-D Fourier transform of the 3-D PSD: 



^■(h+Ah)/ sine p(h+Ah}/ sine 

dO^dOy I dz I dz'K^e~^"'^°'^''^+°'^'^^'>C3D{x,y,z' - z) 

J h/ sin e J h/ sin e 



n p p{h-\-Ah) I sin e p\^n-\-L\n ) / sin e p 

^ J J h/ sine J h/ sine J 
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with ai9i = kiXi . Carrying out the all the integrals except for dfc^, we have 



h 



KSme 



dKPMK,ky,K)^ [^i"^(^fc^AVsine) 



Ksme 



sin£/2Ah A/l^ 

dk'^Ps.D{kx,ky, k'J _ ^ K, 

-sine/2Ah Sin 6 



h/ sine' ft,/ sine' y sine \Ksine 



(Al) 



In the last two steps, the stationary phase approximation is used. Eq. |Alj states that the projected power spectrum 
obeys the same 3-D KT power law, A:"^^/'^, and its RMS is proportional to the square root of the path length through 
the turbulent layer, as would be expected. 
With Pan oc gq jjjj impHes that 



P(a) = i3^in(e)-«/3(a^ + a^)-"/« 
Comparing this expression with eq.[13] in LH2000, we can see that; 

Bl = Ahli^ , 



(A2) 



and has units of mK rad 



-5/3 




Fig. A15. — The coordinate system for the angular power spectrum derivation given in the appendix. 



REFERENCES 



ChamberUn, R. A. & Bally, J. 1994, Appl. Opt., 33, 1095 

Grill, B. P., Ade, P. A. R., Artusa, D. R., Bhatia, R. S., Bock, 
J. J., Boscaleri, A., Cardoni, P., Church, S. E., Coble, K., de 
Bernardis, P., de Troia, G., Farese, P., Ganga, K. M., Giacometti, 
M., Haynes, C. V., Hivon, E., Hristov, V. V., lacoangeli. A., 
Jones, W. C, Lange, A. E., Martinis, L., Masi, S., Mason, P. V., 
Mauskopf, P. D., MigUo, L., Montroy, T., Netterfield, C. B., 
Paine, C. G., Pascale, E., Piacentini, F., Polenta, G., Pongetti, 
F., Romeo, G., Ruhl, J. E., Scaramuzzi, F., Sforna, D., & Turner, 
A. D. 2003, ApJS, 148, 527 

Kovac, J. M., Leitch, E. M., Pryke, C, Carlstrom, J. E., Halverson, 
N. W., & Holzapfol, W. L. 2002, Nature, 420, 772 

Kuo, C. L., Ade, P. A. R., Bock, J. J., Cantalupo, C, Daub, 
M. D., Goldstein, J., Holzapfel, W. L., Lango, A. E., Lueker, 
M., Newcomb, M., Peterson, J. B., Ruhl, J., Runyan, M. C, & 
Torbet, E. 2004, ApJ, 600, 32 

Lane, A. P. 1998, in ASP Conf. Ser. 141: Astrophysics From 
Antarctica, 289-+ 

Lay, O. P. 1997, A&AS, 122, 535 



Lay, O. P. & Halverson, N. W. 2000, ApJ, 543, 787 
Marks, R. D. 2002, A&A, 385, 328 

Marks, R. D., Vernin, J., Azouit, M., Manigault, J. F., & CleveUn, 

C. 1999, A&AS, 134, 161 
Radford, S. 2002, in ASP Conf. Ser. 266: Astronomical Site 

Evaluation in the Visible and Radio Range, 148 — h 

Runyan, M. C, Ade, P. A. R., Bhatia, R. S., Bock, J. J., Daub, 
M. D., Goldstein, J. H., Haynes, C. V., Holzapfel, W. L., Kuo, 
C. L., Lange, A. E., Leong, J., Luokor, M., Newcomb, M., 
Peterson, J. B., Reichardt, C, Ruhl, J., Sirbi, G., Torbet, E., 
Tucker, C., Turner, A. D., & Woolsey, D. 2003, ApJS, 149, 265 

Taylor, G. L 1938, Proc. R. Soc. Lond. A, 164, 476 

Travouillon, T., Ashley, M. C. B., Burton, M. G., Storey, J. W. V., 
&; Loewenstein, R. F. 2003, A&A, 400, 1163 

WoUack, E. J., Devlin, M. J., Jarosik, N., Netterfield, C. B., Page, 
L., & Wilkinson, D. 1997, ApJ, 476, 440 



